Do supernovae favor tachyonic Big Brake instead de Sitter ? 
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Abstract. We investigate whether a tachyonic scalar field, encompassing both dark energy and dark matter-like features will 
drive our universe towards a Big Brake singularity or a de Sitter expansion. In doing this it is crucial to establish the parameter 
domain of the model, which is compatible with type la supernovae data. We find the 1 o~ contours and evolve the tachyonic 
sytem into the future. We conclude, that both future evolutions are allowed by observations, Big Brake becoming increasingly 
likely with the increase of the positive model parameter k. 
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INTRODUCTION 

With the discovery of cosmic acceleration {H the quest for modeling dark energy |01 has started. Besides the most 
simple cosmological constant, other models based on various perfect fluids with negative pressure, like Chaplygin 
gas fl, minimally and non-minimally coupled scalar fields and fields having non-standard kinetic terms were 
advanced. The latter ones include as a subclass the models based on different forms of the Born-Infeld-type action, 
which is often associated with the tachyons arising in the context of string theory |g], Due to the non-linearity of the 
dependence of the tachyon Lagrangians on the kinetic term of the tachyon field, the dynamics of the corresponding 
cosmological models appears to be very rich. 

The tachyon model studied in paper [7] contains a 2-fluid analogue scalar field T ', the dynamics of which is given 
by a simple potential, depending on two parameters, A and k. The model is homogeneous and isotropic. A phase 
space diagram in the tachyonic field and its derivative s = f shows 5 type of distinct cosmological evolutions possibly 
occurring for the model, some of them containing regimes where s is superluminal. All evolutions originate from one 
of the Big Bangs of the model, but they either end in a de Sitter infinite exponential expansion, as the ACDM model 
does, or in a future singularity characterized by a regular scale factor a, vanishing Hubble parameter H and energy 
density e, but infinite s and pressure p. Most notably, the second time derivative of the scale factor goes to — °°, the 
reason why we call this singularity a Big Brake. 

A kinematical analysis |8|] predicted the existence of such singularities, named sudden future singularities. From 
a combined kinematical and observational reasoning alone, sudden future singularities could occur as early as in ten 
million years H, however no underlying dynamics is known to support this. 

Classically the Big Brake singularity is stable. This can be seen by a series expansion of the scale factor in the 
vicinity of the singularity and checking the stability conditions advanced in Ref. [10]. Its quantum study indicated 
singularity avoidance 01 111 . 

Recently 11211 the compatibility of the model with type la supernovae observation has been investigated. After we 
present some basic features of the model in Section 2, in Section 3 we give more details on this compatibility check, 
in terms of the original variables employed in Ref. 10]. Then in Section 4 we stress the crucial difference between 
negative and positive values of the model parameter k. While for the former all evolutions end in the de Sitter attractor, 
for positive k the 1 a contour compatible with type la supernovae contains both states which evolve into de Sitter or 
into a Big Brake. In this dynamical model the Big Brake can occur no earlier than 10 8 years. 



The Big Brake singularity belongs to the class of soft cosmological singularities, which also includes other rep- 
resentants 11311 , Other types of singularities arising in the study of various dark energy models include the Big Rip 
singularity [ 14], present in some models with phantom dark energy [15]. The possibility of existence of a phase of 
contraction of the universe, ending up in the standard Big Crunch cosmological singularity was also considered [16]. 

Unit convention: the Newtonian constant is normalized as EtzG/3 = 1 and we take c = 1. 



THE TACHYONIC MODEL 

We consider the fiat Friedmann universe ds 2 — dt 2 — a 2 (f)t/Z 2 ,where dl is the spatial distance and a the scale factor, 
containing a tachyon field T evolving according to the Lagrangian 



L=-V{T)^\~g^T^T v . (1) 
The energy density and pressure of the tachyon field for the Friedmann background are: 



e = 



VT^t 2 

We shall consider the model with the tachyonic potential 10] 



V(1] p = -V{T)Vl-t 2 . (2) 



V(T) = A r^J 1 - (1 + A:)cos 2 (I ^A(l +k)r) , 

sin 2 ( U/MT+ET ) V V 2 / 



(3) 



where A is a positive constant and — 1 < k < I. The dynamics of the tachyonic field is encompassed in the system: 

t=s, (4) 

* = -3VV(l-s 2 ) 3 / 4 s-(l-. 2 )^, (5) 
while gravitational dynamics is given by the Friedmann equation 

H 2 = e, (6) 

where the Hubble variable H is defined as H = d/a. 

For a negative parameter k, the evolution of the system (0])-(0 is always characterized by — 1 < s < I. The evolutions 
start from a Big Bang and the system has an attractive node at 

T = JL—= , *o = 0, (7) 
3^A(l+fc) 

which corresponds to a de Sitter expansion with Hubble parameter Ho = \/A- (For more details see Ref. id.) 

The case k > is much more richer (see Fig. [T]>. The dynamical system (|4]i-(|5]i has three fixed points: the node (O 
and the two saddle points with coordinates 



T\ = — . 2 arccosW \ — \, s\ = 0, (8) 

3y/(l+k)A V 1+*' 



and, respectively, 



3y/(l+k)A \ ~ V l+k 



T 2 = — ■= I n - arccosA / | , s 2 = 0, (9) 



which give rise to an unstable de Sitter regime with Hubble parameter H\ = y (1 + k)A/2\/k > H( 



FIGURE 1. Phase portrait evolution for k > (k = 0.44). 



The most striking feature of the model under consideration with k > consists in the fact that now the cosmological 
trajectories are not confined to the rectangle given by — 1 < s < 1 and Tt, <T < T4, where 

T3 = — arccos— -, , (10) 
3y/(l +k)A VT+k 

2 ( 1 \ 

Ta = 7T — arrr.ns — (11) 

3y/(l+k)A V y/T+kJ 
are the limits of the domain for which the potential V is well-defined. Indeed, the curvature scalar 

3V(T) (4 -3s 2 ) 

R= v , v — '- (12) 

indicates curvature singularities at s = ±1 except when V (T) — at the same time. This happens at the points P, Q, Q' 
and P' , where as the analysis of [7] shows, there is no singularity and the trajectories can be continuated. In doing so, 
the potential should be redefined by multiplying with i, so that it becomes real, W(T) — iV (T). Then in the energy 
density and pressure yl — s 2 will absorb this i, so that in the superluminal regimes we have 

£= W(r) p = W (T)V^l, (13) 

V* — 1 

both positive. 

In what follows, we briefly display all possible classes of cosmological evolutions existing in the tachyonic model 
with k > 0. First of all, note that in the phase space the reflections with respect to the node point T = Tq,s = leave 
the cosmological evolutions invariant. Thus, it makes sense to study only half of the possible initial conditions in the 
rectangle. This rectangle in the phase space (T,s) should be complemented by four infinite stripes (see Fig.l). The 
left upper stripe (the right lower stripe) corresponds to the initial stages of the cosmological evolution, while the right 
upper stripe (and the left lower stripe) corresponds to the final stages. There are five classes of qualitatively different 
cosmological trajectories. In characterizing them, we shall consider only half of the possible initial conditions taking 
into account the reflection symmetry mentioned above. 

The trajectories of class IV begin in the left upper stripe in the point with coordinates T = 0,s = J i^, which 
corresponds to the singularity of the standard Big Bang type. These trajectories climb to some maximal value of 
s, then turn down and cross the point P, entering the rectangle. Then they leave the rectangle through the point Q' 
entering the left lower stripe. Here, after a finite period of time the universe encounters a special type of cosmological 
singularity, which we call Big Brake. At this singularity, the tachyon field has some finite value, its velocity s tends to 
—00, the cosmological radius has a finite value, its first time derivative is equal to zero, while its second time derivative 



FIGURE 2. Time evolution of the Hubble parameter H(t). 



tends to — °°. The trajectories of class I also begin in the point with coordinates T — Q,s — y ip, however, after 
entering the rectangle they end their evolution in the de Sitter node. They are separated from the trajectories of class 
IV by the separatrix T, which inside the rectangle connects the corner P with the left saddle point. The trajectories of 

class II are separated from those of class I by the curve a, which begins in the point with coordinates T = 0, s = J 
passes through the corner P and ends in the de Sitter node. These trajectories begin at the singularity s = 1, T = Tj n , 
where < Tj„ < T*, T* > T3 and end in the de Sitter node. The separatrix ending in the right saddle, separates 
the trajectories of class II from those of class III. The latter, beginning at s = 1 and T = 7}„, where < Tj„ < T4, 
after crossing the corner Q encounter their Big Brake singularity in the upper right infinite stripe. These trajectories 
are separated from those of class IV by the curve %, which passes through the right saddle point and the corner Q. 
Finally, the trajectories of the last class V begin at s = 1 , T = Tj„ > T4 and end in the Big Brake singularity. The time 
dependence of the Hubble parameter for these five classes is represented in Fig. 2. 

We conclude this section by giving some additional formulae characterizing the different types of singularities 
present in the cosmological model under consideration. In the vicinity of the singularity which takes place at the 
horizontal sides of the rectangle (say, at s = 1), we have the following dependence of the function s on I jj: 

s=l-C(T in )(T-T in )\ (14) 

where 



81 A 2 (j-(l+£)cos 2 ^p^ 
C(Tin)= 32 ^4 3^0^ (15) 



Hence the energy density is 



V{T in ) 



V2C(f~(T -Tin) 1 9(T-T in f 
In the vicinity of the singularity T — T m =t and 



(16) 



9f 2 

while the Hubble variable is 



4 

£ = ^9 (I 7 ) 



3f' 

just like in the dust-filled universe born in the vicinity of the Big Bang singularity. 
For the universe born in the point s = J , T = the potential W behaves as 



(18) 



W ^=wh^ (19) 



where T = J ^~t. The energy density behaves as 



(20) 



9(l+/fc) 2 f 2 ' 
while the Hubble variable is 

H= , 2k r ■ (21) 

3(1+*)/ 

Thus, one can note that the universe has at this point a Big Bang singularity and behaves in such a way as if it were 
filled with a perfect barotropic fluid with equation of state parameter w = 4 . 

We can also describe the behavior of the universe in the vicinity of the final Big Brake singularity following the 
logic of paper [7]. Consider the universe which is approaching the Big Brake in the lower left stripe at some value of 
the tachyon field Tbb- Correspondingly, the variable s approaches — °°. Analyzing Eq. (0 in this limit we have 

' s l = (i«)" 3< '™-' r2/! ' <22) 

where f#g means the moment of Big Brake. Now, using the formula © and the energy density from $13[ , one easily 
finds 

H= (9W^V<\, BB _,y„ (23) 



Thus, we see that when t — ► tsB, the Hubble variable H vanishes while its time derivative diverges, tending to — °°. It 
is important to emphasize that the value Tbb is rigorously positive Tbb > 10]. 



CONFRONTATION WITH TYPE IA SUPERNOVAE 



Following Ref. [18], in Ref. 111211 we have presented in detail how to perform a j 2 -test for comparing the prediction of 
the model with the available type la supernovae taken from Ref. IU7I1 . In order to do this, we introduce more suitable 
dimensionless variables 

H=— ,V = ^r,Q. A = -^,f=H T, (24) 
Ho Hi 

where Ho is the present value of the Hubble parameter. In general, for any variable / (z) we will denote by /o = f(z = 
0). As a follow-up, we also introduce a new tachyonic variable 

y = cos Qv"a(1 (25) 
and switch from the time derivative to the derivative with respect to the redshift z by 

± = -H{l+z)l (26) 
at az 

Then we rewrite the equations ©, ©, © in terms of the new variables H, s, y and perform the ^ 2 -test. For this we 
employ 

d I dt \ 1 



dz V l+zy h ' 

where = H§dh and dL is the luminosity distance for a flat Friedmann universe: 



(27) 



A «-( 1+ «>jfs^) (28) 



The results are represented on the figure panel[3] 



k=-0.4 



k=0.4 




s 




I 



FIGURE 3. The fit of the luminosity distance vs. redshift for k = —0.4 (left) and 0.4 (right). The white areas represent forbidden 
parameter regions where the tachyonic field would be superluminal today. The contours refer to the 68.3% (la) and 95.4% (2o") 
confidence levels. For increasing values of |Jt| < 1 the well-fitting regions are increasingly smaller 1121 . The colour code for x 2 is 
indicated on the vertical stripes. The model is symmetric under the simultaneous change of signs yo — > — yo and .so — > —So, thus 
there is a double coverage of the parameter space. 



FUTURE EVOLUTION 



±°° to finite parameter 



(29) 



In order to avoid the double coverage of the parameter space; also to bring the Big Brake at s 
distance, we introduce the new variable 

1 

w = ~ . 

1 +s 2 

We do this by numerical integration of equations of motion from z = towards negative values of z. We represent the 
future evolution for k = ±0.4 on Fig [4] The evolution curves start from the allowed region (wo,yo) m me plane z = 0. 
The final de Sitter state is characterized by the point (w^s = ^,yds — ®>ZdS = — 1)« me Big Brake final state by points 
(w B b = 0,-1 <y BB <0,-l <zbb<0). 

Whereas all trajectories with k = —0.4 end up eventually into the de Sitter state, those with k = 0.4 can either evolve 
into the de Sitter state or into the Big Brake state, depending on the particular initial condition (wQ,yo). These are 
generic features holding for negative and positive values of k, respectively. In Ref. 1 12] we have also found, that future 
evolutions towards the Big Brake singularity of the universes selected by the comparison with supernovae data become 
more frequent with increasing (positive) k. 



k=-0.4 



k=0.4 





FIGURE 4. The future evolution of those universes, which are in a 68.3% confidence level fit with the supernova data. The 1(7 
contours (black lines in the z = plane) are from Fig [5] (the parameter plane (yo,wo) is the z = plane here). The figures are for 
k = —0.4 (left) and k = 0.4 (right). The coordinate w is related to the tachyonic speed as w = l/(l+i 2 )).Inthe/t = 0.4 case from 
the Ict parameter range the universe evolves either into a de Sitter regime or towards the Big Brake singularity. 



TABLE 1. Properties of the tachyonic universes with k = 0.4 which 

(a) are within 1 o" confidence level fit with the type la supernova data and 

(b) evolve into a Big Brake singularity. Columns (1) and (2) represent 
a grid of values of the allowed model parameters. Columns (3) and (4): 
the redshift z* and time f„ at the future tachyonic crossing (when s = 1 
and the pressure becomes positive). Columns (5) and (6): the redshift 
zbb and time tsB necessary to reach the Big Brake. The former indicates 
the relative size of the universe when it encounters the Big Brake. (The 
values of f„ and tBB were computed with the Hubble parameter Hq = 73 
km/s/Mpc.) 
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For all future evolutions encountering a Big Brake singularity we have computed the actual time f#g it will take 
to reach the singularity, measured from the present moment z = 0, using the equation (H t)' = ~H l (1 +z) _1 ). 
The results are shown in Table Q] The parameter values at which the pressure turns from negative to positive (at the 
superluminal crossing) are also displayed. 

In Ref. [ 12] we have also shown that the Big Brake final fate becomes increasingly likely with the increase of the 
positive model parameter k. 
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